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Abstract. We study the magnetic properties around the Mott transition in the 
Kagome lattice Hubbard model by the cellular dynamical mean field theory combined 
with quantum Monte Carlo simulations. By investigating the q-dependence of the 
susceptibility, we find a dramatic change of the dominant spin fluctuations around the 
Mott transition. The spin fluctuations in the insulating phase favor down to the lowest 
temperature a spatial spin configuration in which antiferromagnetic correlations are 
strong only in one chain direction but almost vanishing in the others. 



PACS numbers: 71.30. -f-h 71.10.Fd 71.27.-|-a 

Geometrical frustration is one of long-standing problems in spin systems. Recently, 
the frustration effects have attracted much attention also in itinerant electron systems. 
The observation of heavy fermion behavior in LiV204 pQ, which has the pyrochlore 
lattice structure with a corner- sharing network of tetrahedra, has activated theoretical 
studies of electron correlations with geometrical frustration |21 El E] ■ The discovery of 
superconductivity in the triangular lattice compound NaxCo02 • yH20 |3] and the 13- 
pyrochlore osmate KOS2O6 [HI has further stimulated intensive studies of frustrated 
electron systems. Geometrical frustration has uncovered new aspects of the Mott 
metal-insulator transition. Among others, a novel quantum liquid ground state was 
suggested for the insulating phase of the triangular lattice , and this may be relevant 
for frustrated organic materials such as k-(ET)2Cu2(CN)3 p]. 

The Kagome lattice (Fig. H} is another prototype of frustrated systems showing 
many essential properties with the pyrochlore lattice. It is suggested that a correlated 
electron system on the Kagome lattice can be an effective model of NaxCo02 • yH20 
by properly considering anisotropic hopping matrix elements in the cobalt 2>d orbitals 
IHj. The electron system on the Kagome lattice in the metallic regime was studied 
recently by using the fluctuation exchange (FLEX) approximation |T^ and quantum 
Monte Carlo (QMC) method [TT], etc OCSI. In our recent paper [T3j, we have studied 
electron correlations in the Kagome lattice Hubbard model, and found the first-order 
Mott transition at the Hubbard interaction U /W ~ 1.37 {W: band width). 
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Figure 1. (a) Sketch of the Kagome lattice and (b) the effective cluster model using 
three-site cluster CDMFT. (c) First Brillouin zone of the Kagome lattice. 



In this paper, we focus on the magnetic properties of the Kagome lattice Hubbard 
model. By applying the cellular dynamical mean field theory (CDFMT) ^3], we 
discuss the effects of geometrical frustration around the metal-insulator transition where 
frustration is stronger than in the weak coupling regime. We consider the standard 
Hubbard model with nearest-neighbor hopping t > on the Kagome lattice, 



u 



(1) 



with rijo- = c-g-Cjo-, where c]^ (cj^j) creates (annihilates) an electron with spin a at the 
site i. In the following, we use the band width W = 6t as the energy unit. The 
dynamical mean field theory (DMFT) jTHj has given substantial theoretical progress 
in the field of the Mott transition but it does not incorporate spatially extended 
correlations. Therefore in order to take account of geometrical frustration, we use 
CDMFT, a cluster extension of DMFT [13 El CH] , which has been successfully applied 
to frustrated systems such as the Hubbard model on the triangular lattice [ini EOl 121] • 
In CDMFT, the original lattice is regarded as a superlattice consisting of clusters, 
which is then mapped onto an effective cluster model via a standard DMFT procedure. 
As shown in Fig. ^ the Kagome lattice Hubbard model is mapped onto a three-site 
cluster coupled to the self-consistently determined medium, 
r-/3 

J2 4- (^) ^isi 



^eff = / drdr' J2 (r) G't (r - r') c^. (r') + U j drj^n,^ (r) n,^ (r) . (2) 







Given the Green's function for the effective medium, Q„, we can compute the cluster 
Green's function G„ and the cluster self-energy So- by solving the effective cluster model 
with QMC method [22] ■ Here, Q^j, G^, and Sg. are described by 3 x 3 matrices. The 
effective medium is then computed via the Dyson equation. 
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where /i is the chemical potential and i (k) is the Fourier-transformed hopping matrix 
for the superlattice. 
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Here the summation of k is taken over the reduced Brillouin zone of the superlattice (see 
Fig. ^c)). After twenty iterations of this procedure, numerical convergence is reached. 
In each iteration, we typically use 10^ QMC sweeps and Trotter time slices L = 2W/T 
to reach sufficient computational accuracy. Furthermore, we exploit an interpolation 
scheme based on a high-frequency expansion of the discrete imaginary-time Green's 
function obtained by QMC in order to reduce time slice errors. 

We now investigate the magnetic correlation around the Mott transition in the 
Kagome lattice Hubbard model. We calculate the wavevector-dependence of the static 
susceptibility. 



/" / \ 

X^5(q) = / drY^ ^4^^ (r) Ck+q^^ (r) cl^^g^ (0) c^'sr (0)^ , 

"^0 k,k' 



(6) 



where 7,5 = 1,2,3 denote the site indices in the unit cell. We employ the 
standard procedure in DMFT to calculate x-rsi'i) [IB], which includes nearest-neighbor 
correlations as well as on-site correlations. In order to obtain x-ys{^), we first calculate 
the two-particle Green's function in the effective cluster model 

rP rP rP rP 

C^S {i^h i^m) = ^ / / / / dTidT2dT'idTi 

Jo Jo Jo Jo 

(7) 
(8) 



Jo JO Jo Jo 

C^5{ri,T2,T3,n) = (T^c\^{Ti)c^i{T2)cl^{T3)cs^{n) ^ 

and extract the vertex function F^^ {iui, ium) via the Bethe-Salpeter equation. 



F = CO 



-1 



c-\ 



where is the bare two-particle Green's function. 
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On the other hand, the bare q-dependent Green's function in the lattice system is 
calculated by 



C?5(q : i^i) = -fYl ^^^i (k + q : iuJi) g^s^ (k : iui) . 



(11) 



By using Eqs. © and (fTTjl . we can compute the lattice q-dependent Green's function. 
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Taking account of the phase factor, we finally obtain the q-dependent susceptibility. 



x^siq) =t2^C^5 (q: 



(13) 



It is convenient to introduce Xm{^} for three normal modes (m = 1,2,3) by 
diagonalizing the 3x3 matrix X7<5(q)- the upper panels of Fig. |2l we show the 
three eigenmodes of the susceptibility at T/W = 1/30 for several values of interaction 
strength U/W. In the noninteracting case, the largest eigenvalue of the susceptibility 
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Figure 2. The wavevector dependence of the static susceptibility Xm{(l) for several 
values of U/W at T/W — 1/30. The three dimensional plots of Xm(q) are shown in 
upper panels, from top to bottom to = 1, 2, 3. The two dimensional plots in the lower 
panels show the dominant mode of the susceptibility upper panels. Hexagons 

in figures denotes the first Brillouin zone as shown Fig. ^(c). 




Figure 3. The enhanced spin correlations in the insulating phase U/W = 1.4 at 
T/W = 1/30. 

Xi(q) takes a maximum at six points in the Brillouin zone, the second largest one 
X2(q) has a maximum at q = (0,0), and the smallest one Xi{^) takes maxima at 
the corners of the Brillouin zone. As U/W increases, localized moments are formed 
and the susceptibihty is enhanced as expected. In particular, the q = (0, 0) peak of 
X2(q) becomes strongly enhanced, which is consistent with the previous QMC study 
[TT] . On the other hand, the dominant mode of the susceptibility shows only 

weak q-dependence. As U/W increases, is enhanced not only at the six points 

above mentioned but also enhanced on the lines passing through F and M points, so 
that the q-dependence of the susceptibility gets suppressed and becomes much weaker 
at U/W = 1.1 than in the noninteracting case. This behavior is consistent with the 
previous FLEX calculation in the weak coupling regime ■ We confirm that the feature 
of the suppressed q-dependence of the dominant magnetic mode due to geometrical 
frustration persists up to fairly large-[/ region. 

We further find notable results in the insulating phase. The Mott metal-insulator 
transition occurs at U/W ~ 1.37 fl]. As shown in the lower panels of Fig. El 
once the system enters the insulating phase, the q-dependence of dramatically 
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changes its character due to the enhancement of short range antiferromagnetic (AF) 
correlations [14 . At U /W = 1.4, the susceptibihty takes the maximum value along 
the three lines in q space instead of the six points in the weak coupling regime. 
Furthermore, by investigating the eigenvectors of Xi{^)j we find that two spins in the 
unit cell are antiferromagnetically coupled but the other spin is free. Therefore, at 
these temperatures, the enhanced spin fluctuations favor a spatial spin configuration in 
which one-dimensional (ID) AF-correlated spin chains are independently formed in three 
distinct directions. The three types of enhanced spin correlations are illustrated in Fig. 
01 This spin correlation is one of the naturally expected spin correlation on the Kagome 
lattice, because it stabilizes antiferromagnetic configurations in one direction, which is 
more stable than the naively expected spin configuration having a singlet pair and a 
free spin in each cluster. These ID correlations in the finite-T Mott insulating phase 
are different from the results for the Heisenberg model on the Kagome lattice with the 
nearest-neighbor exchange obtained by both classical and semi-classical approximations 
[23 ; but are similar to the q = structure predicted for the classical Heisenberg 
model with a further neighbor exchange [21]. The essential difference from Ref. is 
that there is almost no correlation between the different chains in our results for the 
Hubbard model. It remains an interesting problem to compare the q-dependence of 
the susceptibility in the large-t/ and low-T regime with the results for the Heisenberg 
model. 

Finally, we show the maximum value of the susceptibility Xmax as a function of 
U/W in Fig. |3] As U/W increases towards the Mott transition point, 1/Xmax decreases 
almost linearly as expected. When U /W is further increased and the system enters 
the insulating phase, the spin correlation changes into the ID one as discussed above, 
and then Xmax saturates. Therefore, we find no evidence of the real instability to an 
ID ordering in the present calculation. However, such enhanced spin fluctuations affect 
low-energy dynamics in the insulating phase ^1] . 
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In summary, we have studied the magnetic properties around the Mott transition 
in the Kagome lattice Hubbard model by means of CDMFT combined with QMC We 
have investigated the q-dependence of the susceptibility, and obtained consistent results 
with the previous studies: the second-largest eigenmode of the susceptibility shows 
the strong q-dependence, taking a maximum at q = (0,0), and the q-dependence of 
the maximum eigenmode of the susceptibility becomes suppressed as U increases in the 
metallic phase. We also find a dramatic change of the dominant spin fluctuations around 
the Mott transition. The spin fluctuations in the insulating phase favor a novel spatial 
spin configuration in which ID AF-correlated spin chains are independently formed in 
three distinct directions. Although the spin liquid state or other nonmagnetic ordered 
states may be stabilized at zero temperature the enhanced ID spin correlations 
could emerge in the finite-T Mott insulating phase. 
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